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* 

Donald P. Gaver 
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1 . Introduction . 

Many situations suggest the study of the following class 

of models: a point stochastic pacing process , P, with counting 

process {N(t),t^0} develops over time, and associated -with each 

event or "point" is a real-valued random variable, X; let X n 

denote the n— such random variable, realized at the moment 

N ( t) first reaches n. Here P may be interpreted as specifying 

the occurrence of demands — or supplies — of goods, water, payments 

or reimbursements, opportunities, e tc . , while X represents a 

generic demand or reward size, X n being the magnitude of the 
th 

n— event. Alternatively P governs the occurrence of system 
shocks, possibly caused by floods or earthquakes, and X is a 
generic shock magnitude. 

The objective of this paper is to introduce various simple 
models for P, and for (X n ) and to study the occurrence of 
record events ("new highs") in the X-sequence as they appear under 
the stimulus of P. Related investigations are those into the 
classical record problem by Barton and Mallows [1] , Chandler [2] , 
Foster and Stuart [7], and Shorrock , e.g. [13], [14]. The studies 

of extremal processes by Dwass [4], [5], see also Resnick [12] 

and Resnick and Rubinovitch [11] also are pertinent. In particular 

* 

Research support from NSF Grant AG476 at the Naval Postgraduate 
School is gratefully acknowledged. I wish to acknowledge useful 
comments by George Humfeld, Mark Westcott, and the referee. 



Pickands, in [10] , allows record events to be Poisson paced. The 
results reported in this paper are mainly explicit forms, and 
domains of attraction, for first record-time distributions when 
P is one of several familiar point processes (Poisson, renewal, 
pure birth) . In addition we point out that our distributions may 
be interpreted as equipment lifetime distributions in the relia- 
bility sense, and that the classical "marriage problem" of sequen- 
tial decision theory, cf. De Groot [3], may be viewed as a record 
problem . 
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2 . Process Specification . 

A point stochastic pacing process , P, with counting 
process {N(t),t^0} governs the occurrence of values , denoted 
by {X , n=0 , 1, 2, . . . } . Assume P and {X } to be mutually 



n 



n 



independent, and the {x^} to be an i.i.d. sequence with abso- 
lutely continuous d.f. F(x). Mainly we shall assume that P 

is an orderly process, so we can speak of the times between 

th 

successive events in P; let the n— of these be (T , n=l,2,...} 

n 

so that the time of the n— process event is t = . T. , 

r n ^ 1=1 l 

n = 1,2,... Suppose that at t = 0 a reference value, Xq , 

is available; Xq is a random variable with d.f. (x) . 

Definition . A first record with respect to X Q occurs at 

t = t if (a) X. < x rt , all i < n, (b) X > X A , in which 
n^ l U 1 n^ U 

case T, = t = the first record time, and X = R. , the first 
In. n. 1 

-L th 

record value. Subsequently, a k— record occurs at time 



= l . . r . if (a) X. < R. . , 

n,. ^ 1=1 l l k-1' 



*k 

X > R 



n k-i ' 1 < v (b) 



, . , and t. = t - t , R, = X Alternatively 

n k k n k n k-l k n k 



x. = inf{t: max X > X A } 

1 , . > n 0 

n^N ( t) 



( 2 . 1 ) 



R 1 X N(t 1 )' 



( 2 . 2 ) 



T.. + T_ = inf{t: max X > R. } 

_L 2 / » \ n -L 

n^N ( t ) 



R 2 X N(t 1 +t 2 )' 
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3 . Poisson Paced First Records, with Variations. 



Suppose P is Poisson with time -dependent rate X (t) , 

and hazard A (t) = [ X(t')dt'. If F n ( * ) denotes the d.f. of 

J 0 

Xq, and F(*) the d.f. of {X n , n=l,2,...}, then 

c z 

P(x (1) £ ' R 1 G = | e ^ F dF^ (x) dA (t) dF (z) (3.1) 

from which various special cases of interest immediately appear. 



3.1 F Q ( • ) = F ( • ) . 

In this case the integration is immediate, and 

-A(t) (l-F(z) ] -A (t) 

P{t 1 e (dt) ,R ± e (dz) } = S 7 * dA(t)dF(z) . (3.2) 

The marginal d.f. of record size is obtained by integrating out t 
in (3.1) ; the density is 



P{R 1 e (dz)} = -log [1-F (z) ]dF (z) 



(3.3) 



just as in the classical setup. The marginal d.f. of x^ is given 
by integrating first on z. Most simply, however, 



P{x 1 > t | X Q = x} 



-A ( t) [1-F (x) ] . 
G / 



(3.4) 



Removal of the condition using 

Ptx-L > t} 



F Q = F yields 



1 - e 



-A (t ) 



A ( t ) 



(3.5) 



Note that in this, and in other, cases P{x^ > t} is given by 
integration of the generating function of the counting process of 
P: 

P{x.>t}= f l e “ A(t) - ( - ;=y F n (x)dF(x) . (3.6) 

J n=0 n * 
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Further note that since A(t) is non-decreasing, then if h > 0 



E[x 1 ] = 



f°° 


rh -A ( t ) r 


P{t. > t }dt = 

J 0 1 J 


i/Vu, dt+ | 



1 - 



,-A ( t) 
Tit) 



dt. 



so 



(1-e 



-A(h) 



dt 

A (t) 



£ E [t x ] £ h + 



r 

•* h 



dt 
A (t) 



E [x 1 ] < 00 iff 

For time homogeneous Poisson P it is clear that no finite moments 
exist . 



dt 

TCtY 



V h > 0. 



(3.7) 



3.2 Random (Gamma) Hazard, f q ( ' ) = F (‘). 

If {A(t)} is a gamma process, then E [e ^ 
and {n ( t ) } is a negative binomial process. Therefore 



P{t > t} 
1 



and 



rn 



y+l-F (x) 



kt 



dF(x) 



y 

l-kt 




1-kt 




(3.8) 



p(t > t } -H_ = 1 

1 1 * kt E [A ( t) ] 



as t 



(3.9) 



Letting k decrease but such that — = 1 induces an increasingly 
long Pareto tail, and a buildup of mass near zero. A similar 
effect is achieved by mixing Poisson P's rate parameter by the 
gamma (applying a gamma prior to X) : 



P{t 1 > t} = 



r<*> -Xt . , . . k-1 

1-e „-y X (yX _ y 

0 Tt e " T (k)" MdX - t 



-k+1 

C 1 + 7r) " 1 



-k+1 



(3.10) 



If P has been observed throughout the time interval (-T,0), 

and j events occurred, then the posterior is gamma with y = T+ a, 

— Xa ( Xa ) k - 

k = j + b, where the prior density is p(X) = e — TTb) — a an< ^ 

(3.6) provides a Bayes-type predictor of the distribution of time 
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until the first record, the latter being referred to an initial 
Xq with d.f. F. 



3.3 Process Observed for Initial Time T. 

If the P-paced process has been observed for (-T,0) and 
further records are referred to the largest X occurring therein, 
then 



F„(x) = e - A < T)tl - F(x)) , 



(3.11) 



and application to (3.4) yields 



P{t ± > t} 



A (T) n _ - [A(T)+A(t) ] , 

A(T)+A(t) 1 



(3.12) 



with confidence given by the r.h.s. of (3.12) there will be no new 
maximum in (0,t) greater than that in (-T,0). 



3.4 The Distribution of t . 

n 

Since t , the time between n-1— and n— records 
n 

does not depend upon the distribution of X, we may capitalize 
upon the Markovian property of the exponential, F(x) = 1 - exp(-x), 
to characterize the distribution of t for any n. Now for the 
time-homogeneous Poisson and exponential X, 

. -x 

P{t e (dt) |R i = x } = e e Xe x dt. (3.13) 

n n-1 



We can write the characteristic function (ch.f.) of log x as 

n 

follows. First, 
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E [e 



iz log t 



n 



R =x] = 
n-1 



• \ a. ~x 

12 -A te . -x , 
t e Ae dx 



i 12 

rl'v izx 

= irJ e 



iz -u , 
u e du . 



(3.14) 



Next remove the condition using the fact that R n _^ is gamma: 



iz log t . 

E[e "] = (£•) 



IZ 



izx -x X 
e e 



n-1 



(n-1) ! 



dx 



iz -u , 
u e du 



_ e iz log A f 1 



n 



(to) 



iz log u -u , / o i 

: ^ e du, (3.15) 



0 



which may, on the basis of ch.f. unicity, be interpreted as saying 
that the distribution of Ax n is that of the sum of two independ- 
ent random variables, one being gamma with mean and variance 
equalling n, and the other .being the logarithm of a unit expo- 
nential. Moreover, the ch.f. of — log At is 

n ^ n 



log At 



iz 



n 



E [e n ] = (1-^) 



-n 



iz -log u 

n -u , 

e e du 



(3.16) 



e" 



iz 



as 



n 00 



and invocation of the continuity theorem for the ch.f. shows that 

n /TF -> e in probability. Finally, consider the ch.f. of 
log At - n 

. As n becomes large the gamma component dominates 

/n 

and the continuity theorem shows that the latter sequence of r.v. 
tends to N(0,1) in law. These results are analogous to those 
of Neuts [9] for the classical record problem. 
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4. Renewal-Paced Records. 



Suppose the pacing process P is an arbitrary renewal 
process with i.i.d. interevent times (U n ) having d.f. H(*)» 
and H(0+) = 0. Then 



P{t 1 > t | X 0 = x) = l [H n *(t) - H (n+1) *(t) ]F n (x) (4.1) 

n=0 



and, Laplace transforming, the following explicit form appears: 

e -st P{i, > t | X. = x}dt = 1 -- ^ (S) ^ (4.2) 

J 0 s 1-H (s) F (x) 

A 

where H is the Laplace-Stielt jes transform of H. If F Q = F, 
we integrate to obtain 



00 

e _St P{t 1 > t}dt = 1 - ~ I ? (s) log [1-H (s) ] , (4.3) 

•'0 -sH (s) 



SO 



-ST 

E [e 



1 + log [1-H (s) ] . 

H (s) 



(4.4) 



Example 1 ♦ If P 
P{T 1 >U ~ Xt as 



is time homogeneous Poisson we have seen that 

• • . . . (lforts 1/X 

just as is true when H (t) = < 



t + 



0 for t < 1/X, 



i.e. when records are regularly paced. 



Example 2 . If H is a stable distribution, then H(s) = e 

(0 < a < 1) and this and (4.3) yield 

a 

r°° . s a 

e _st p{t > t> = - — — log [l-e -s ] - -as a log s. (4.5) 
J 0 1 s 



so by a Tauberian theorem (Feller [6], p. 447), 

Pfi^ > t} ~ r ' J 1 ' - ' a )~ a;Lo 9 t (t->-<») (4.6) 
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1 - 1/2 
For instance if a ~ 2 then H tails off like t while 

the d.f. of tails off like t log t. 

4.1 Domains of Attraction. 

Suppose that n independent copies of accumulated 

n 

and summed: S = £ T.. (i) . We record several facts about the 

n i=l 1 

behavior of S n for large n ascertainable from general theory; 
see Feller, II [6], and Gnedenko and Kolmogorov [8], p. 175, 
Theorem 2 . 

A. If P^ T a > t) ~ as is true when P is Poisson, 

regular-spaced, and in other situations as well, S n is attracted 
to a stable law of order a = 1. Direct expansion of the rele- 
vant Laplace transform around s = 0 shows that s n / n log n -*■ 1 

in probability. 

B. If P is renewal with stable law of order a (0<a<l) 
interevent times, then because of the asymptotic behavior of 4.6), 
the further results of Feller [6] , pp. 448-449 immediately imply 
that 



— ( $ ) z (4.7) 

a a 

n 

Z a being once again stable of order a . Here a^ satisfies 

a a ~ an log a (4.8) 

n ^ n 

1/a . 

and a n > n for large n: obviously a strong norming is 

required to bring the stable-paced record times back to stable. 

A suitable normalizing constant derived from (4.8) is seen to be 
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*- (n log n) 



1/a. 



(4.9) 



This result also follows immediately from a transform continuity 



theorem. Simply scale by a n , (a n +») . 



S x,-n n 



f = l ]) 



= (1 + [exp(^-) a -l] log[l - exp(-(^-) a ] } n (4.10) 

cl a. 

n n 



and expand in Taylor's series to see that if 

a 



n log a 



the above (4.10) tends to e 



-s 



a 

i 

n 



-* 1 then 



the stable law (a) transform. 
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5 . Furry Records . 



For variety let P be the simple pure-birth or Furry-Yule 
process ([6], I, p. 450) generated by X n = nX. Start with a single 
progenitor at t = 0 (of value Xq drawn from F(x)). Now N(t), 
the number of new births in (0,t), is geometric, and each is 
endowed with a value, X, that is an independent copy of X^. 

Hence 

N ( t) e""^ 

g(z,t) = E[z Wir; ] = ^ rr- (5.1) 

1-z (1-e ) 



and 



P{T 1 >t) 



| g (z, t) dz = 
0 

-Xt 



Xte 



= Xt 

" -Xt At 7 

l-e e -1 



e Xt dF(x) 
l-F(x) [l-e" Xt ] 



00 



= l 

k=0 



B. 

l 

rr 



(At) 1 , 



(5.2) 



where the are the Bernoulli numbers. Integration of (5.2) 

yields the mean of the first record time: 



Elr,) = X’hl + ^ + jr +...]= i ^ « 3^65 . (5.3) 

Thus the Furry P process encourages a record to occur after a 
mean delay of very nearly half way between the mean second and 
third Furry jump times. 

5.1 General Pure-Birth P; Laplace Transform. 

If P is pure-birth with general intensity X n we find 
by following a recurrent events argument and an integration that 
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00 



00 



( 5 . 4 ) 



J 

0 



e" st P{T 1 >t}dt 



n=l 




TT 

j=i 



x . 

J_ 

X . +s * 
J 



Now apply a Tauberian theorem allowing s -*■ 0 to obtain 



E[t 1 ] 



I 



n=l 




/ 



which shows that the first record time possesses a mean iff 
series on the r.h.s. converges. 



( 5 . 5 ) 

the 
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6. Models. 



The P-paced record process suggests several specific 

models . 



6.1 A Shock Model in Reliability. 

Suppose X n represents the magnitude of the n— shock 

delivered by P. Let Xg represent the "strength" of a component 

manufactured to inhabit the P environment; the component survives 

to time t if no shock in time t exceeds X Q . Hence (3.5), 

(3.8), (3.10), and (3.12) all may be interpreted as component life 

distributions. It may be reasonable to replace F by 

X 0 

k+1 u 

[F(x)] , k > 0, arguing that the manufacturer will attempt to 

build in a safety factor; safety increases with increasing k. 

In this case (3.2) becomes 



P{T 1 >t} 



(k+1) 



-A(t)w,. ,k. 

e (1-w) dw 



( 6 . 1 ) 



Example ♦ If F (x) is taken to be exponential with 
parameter 0, so E[X] = 0 ^ then, having observed an X we 
can consider setting Xg = r ^X, r < 1 representing a safety 
factor. Then 

_ £3 y 

p{x 0 ^x> = 1 - e (6.2) 



and 



P{x 1 >t} 



-V . r “0X, Q 

-At[e ] „ -r0x, 
e J r0e dx 



r(r+l) 

(At) r 



At 



0 



e 



z 



r- 

z 



1 



r (z) 



dz 



r (r+l) 
( At ) r 



(6.3) 
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as t -*■ 00 , and an exceedingly long tail results if r < 1, 
representing conservative design. The remarks of Section 4 
indicate that sums of times to failure are attracted to the 
stable law of order r if 0 < r < 1. 

6.2 The Marriage Problem. 

At the moments of event occurrence in a Poisson process a 
decision maker (art collector, suitor, employer, or what have you) 
is shown objects of value X. Allowed to inspect for only a 
finite time, and thwarted from any return to previous opportuni- 
ties, he seeks to select that X of maximum value to occur in 
the interval. 

Suppose the decision maker sets aside a time T to look 
over the field, leaving himself U for choice ; and enjoyment or 
regret. U + T is fixed. Let his decision rule be to examine 
all X-values to occur in T, and then to select the first X- 
record, if any, thereafter. He wishes to choose U in such a 
way as to maximize his probability of ending up with the greatest 
X-value to occur during U + T. Following this rule enables him 
to select the best X to occur throughout (0,U+T) with proba- 
bility nearly e "*" = 0.368..., exactly as in the classical formu- 
lation, cf. De Groot [3], in the event that opportunities come 
thick and fast, i.e. as X -*• <». For suppose he confronts the 
time period (T,T+U) , committed to picking the first record 
therein, if any occurs. Given x was the record during (0,T), 
the probability that the first record to occur subsequently during 
(T,U+T) is also the last is 
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<p(U,x) = 



U 00 

e -Xt[l-F(x)] xdt 



X 



dF(y) e - (U - t)[ 1 - F( y»l 



u 



e -Xt[l-F(x) ] _ e -XU[l-F(x) ] 
U - t 



dt . 



The best opportunity during (0,T) has distribution 



F q (x,T) = e 



-XT [1-F (x) ] 



(6.4) 



(6.5) 



and so the probability that the best overall is picked is 



<Pi (T) = | <p (U,x) dF Q (x,T) 



( 6 . 6 ) 



Now interchange integration order and invoke (3.12) to obtain 



U 



n (T) = J isk { 



T+t 



(l-e- X(T+t) ) 



T+U 



(1-e 



-X (T+U) 



)} (6.7) 



If X ->■ °o the integral tends (bounded convergence) to 



U 



<P(T) = 



T+U 



dt _ T , T+U 

log 



T+t T+U 



= -T log T if T + U = 1, (6.8) 

and differentiation shows that <p is maximized for T = e 1 . 
Finally, <p (e ^) = e agreeing with the classical solution, 

see De Groot [3] , p. 331. Now if the decision maker values the 
time to enjoy his prize as well as its attainment then he may be 
tempted to optimize max^-r^O), T ^ being the time of the first 
and only record to occur after T, measured from T as origin. 
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Then his expected enjoyment is 



<MT) 



U 

T 

T+U 

0 



(U-t) 

T+t 



dt = T log 




UT 

T+U 



= T 




UT 

T+U 



= -T log T - T(l-T) (6.9) 

with maximizing value close to T = 0.2, and i/>(0.2) = 0.16. If 
our decision maker chooses the latter strategy he wins with approx- 
imate probability 0.32 rather than the 0.37 probability of success 
achieved under the classical decision rule. On the other hand, 
if he chooses to observe for T = e ^ his payoff under the time- 
weighted utility is 0.14. Perhaps T = 0.25 will appeal as an 
easily remembered compromise that nearly maximizes both goals men- 
tioned. Finally, if X becomes small it may be seen that optimum 
T -*• 0.5 if the above rule is followed. However, this rule cannot 
be optimum under such light traffic situations: for very small 

X one does best by choosing the first opportunity to appear. 
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